contribution of the diffuse scattering to the measured integrated re- 
able values (up to 30 % for the highest reflections), and thus cannot 
been thought until now. 


Risesies reflection can, as to the dependence of the Bragg angle v9, be written 
in2 


tegra “a ‘Bragg reflection times a factor exp (const = ), due to the effect of the 


Pe 


2 
‘This fact explains why good agreement can be obtained between theoretically 
values and those obtained from intensity measurements and experimentally deter- 
‘bye-Waller factors. | 


Introduction 


The effect, of the thermal motion on the scattering of X-rays by crystals has been 
é igated by Debye [3], Faxén [4] and Waller [5]. In his dissertation Waller gives 
- complete formula for the scattering by small crystals. In a later work [6], he has 
‘ived an expression for the integrated intensity of the diffuse maxima, coinciding 


4 with the Bragg maxima. The present calculations are to be considered as an extension 
of Waller’s paper. 
j _ For a simple cubic crystal, containing a large number of unit cells, but yet so 
a small that the absorption within it may be neglected, the following formulas for the 
__ intensity are valid: 
E all a J, ae Ape 

; I be 

J, =F lFP a: (1) 
4 ; I 2,-2M % My 
4 If 2 — R | F | é 2 ye 
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Fig. 1. of 


P wy 


I is the intensity of the incident monochromatic, plane wave, J the scattered inter 
sity. J, is the Bragg intensity modified by the temperature factor e~*™. J, represe 
the diffuse temperature scattering. R is the distance from the crystal to the poin 
where the intensity is observed. F is the atomic scattering factor. T, is the Laue 
interference function. 7’, represents scattering where only one elastic wave is in- 
volved. It has pronounced maxima at the reciprocal lattice points. The scattering 


corresponding to 5 T, forms a nearly uniform background with only slight maxima 
v=2 & 2 


and may thus be neglected for our purposes. Ts 
In Fig. 1 the sphere of reflection is drawn for a direction of incidence PO. The 
lattice wave that gives the contribution represented by T, to the scattering in the 


direction PQ has the wave vector 20 MQ = 22g, M being the reciprocal lattice point 
lying nearest to Q (Faxén-Waller vector condition). 

Our task is to calculate the energy, corresponding to J,, received by a slit, placed in 
the direction of Bragg reflection, when the crystal is rotated through the Bragg 
position with uniform angular velocity, around an axis, lying in the reflecting crystal 
plane and perpendicular to the wave normal of the incident wave. From the formulas 


(72), (65) and (62”) in Waller’s dissertation we get the following expression for 7',, 

valid for small values of 0: 3 
kT N A 

=> ————_ — > 2) 

1 a Oo B ( ) 


where ; 
A = (% = 0%)? [By (Coy + yg) B27? + Cag (Cry — Cay)(B2 + 72) + 024] + cyel. perm. — 
— 2(% = a 9)(B — Bo)(Cre + C44) &B (by 7? + C44) + cyel. perm. 
B= bib, aR 72 + by (Cy, + Cyy) Ogg (&2B? + B27 + Pa?) + C1 Che | 
4 =the wave length of the incident X-rays ; | 
x =27/A | 


k = Boltzmann’s constant 
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T= temperature 


ON = number of atoms in the crystal 

a =distance between adjacent atoms in the lattice 

0 =22a/i =2n2a|g|, where / is the wave length of the elastic wave 
«a, B, 7 are the direction cosines of g referred to the crystallographic axes. 


ee) Bor Vos and «, f, y are the direction cosines of respectively the incident and the 
scattered X-rays referred to the same axes. 


7 Dy = Cyy — Cyg — 2Cqq, 5g = Cyy + 2CqQ + Cay. 


q C11, Cya, ANG Cy, are the elastic constants. 


Fig. 2 shows a section of the region in the reciprocal lattice space, through which 
- the integration of 7, has to be performed. The direction of the incident X-rays is 
- changed from P’O to P’’0 during the rotation of the crystal. Thus the sphere of re- 
flection moves from AD to BC nearly parallel to itself in the neighbourhood of 
the reciprocal lattice point M. The angles A P’D and BP’ correspond to the angular 
width of the slit. The height of the slit corresponds to an extension of the region of 
integration perpendicular to the plane of the section in Fig. 2. Thus the endpoint 
_ of the vector g moves in an approximately parallelepipedal space, when the crystal 
_ is rotated. To perform the exact integration of 7’, would be a very lengthy numerical 
work. The actual calculations are made according to the following approximate 


- method. 
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yet rH Fi pe eu (2 | 


Introducing polar coordinates & = sin 0 cos $, B = sin 0 sin 4, 7 =cos 0 and. 
(=e)? a (9 — Pp) ae Yo)” i we Sas 


fe 


When the integration with a to ¢ has been performed, the following expression 
for XK is obtained 


A C 
K= —— arct y/o 
Von oh 


(BD-AE+BOx+AC#—-AC#)de 
(Ox? + D) tir Ee OB DD EO eC ae : 
+(2C—5 D*—3CD)2°+ (83CD+3D)x°—-CDe" 


+2 


0 
x =cos0, 
A= by (Cy; + C12), C = by*bo, E = ¢y, Cis, 

B= Cy4(2¢y1 + C4q), D = dy (Cy, + C42) Cae 
A simple formula derived by taking the mean values of each of the terms in the | 


numerator and the denominator of (3) separately gives a good approximation of X. 


a es B Dy (C41 + eye) + Cag (2 C43 + Cay) 
mPPt a5 OT ba + $ By (C11 + C49) Cag + Cy Ch 


7 
Table 1 j 
a ee { 
KC NaCl | 
T° abs 
E 
Kx 10" | Kappr x 10 Bs Kx 10% | Kapp x 104 ae | 
‘oO ‘oO 
A ; 
86 2.228 2.112 —5.2 1.418 1.403 os | 
280 2.464 2.368 —3.9 

290 2.480 2.385 —3.8 1.594 1.588 -0.4 | 

) ee eee 


stal, ol slit ‘and « p is 5 the Eolehe of sts slit, 
ion for , 0, valid for small ne of u, w, and y is given by 


“ ery ’ Seiad 


_ = 2a? (Lu? sin? v9 —4u w sin? vy + w? + y? cos? v). 


er 


surface clanarit of the slit is R? cos vydwd g. 


= Uz, We ; 
=] | | cos vy dudwdy 
oO =S=— —_—_— = — _ 
4u? sin? v)—4uw sin® vy + w+ g? cos? v4 

Uy Wi 
_ and putting eS oe K— a a T. 
4 we get the following expression for the integrated diffuse intensity 
. 8,07 [Fe 
- T 2 KiacbV (6) 


p vaitten in conformity with Darwin’s formula for the integrated Bragg reflection 


rf E,wo |FPe -2M 48 & 

of A? gin 20, 
is the angular velocity of the crystal during the rotation, A is the cell volume =a’, 
and 6V = Na? is the volume of the crystal piece. H, and £, are the total energies 


corresponding respectively to J, and J, entering through the slit during the rotation. 
_ Assuming the slit very long and thus putting 


weineins A ei) nme . 


NY 
A 


; Oe ar Wo Pe, 
 o is given by 


We 2 uy sin” vy 2u, sin” v, 
0 2 sin? % i Wy i We he 
(8) 


w 2u, sin? v 2u, sin? v w 
Fag Ot ee |e (- 2): 
2 sin* vp Wy Wy We 


e/ Here : 
- f(x) =a In [(1—2) sin vy + V2? — 2 sin? v, + sin” CA ete 


+sin v, In [#—sin? vy + V2? — 2a sin? v9 + sin” v1: 


ee 
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obt ; g 6V with (2)~ im (7), w being the ordinary line: 
coefficient of the crystal, and I being the total energy incident per unit . 
crystal. In the comparison below we use the expression (9) for «, which means 
we assume the same orientation for all the elements in the crystal mosaic. 


Measurements made by James and Brindley [1] : 
ie 


Table 2 gives the value of « for the reflections of KCl used by James and Brind ley 
to determine # in the expansion 


pr +242.) (10) 
deduced by Waller [5]. In this expression y does not depend on the atomic forces, 
and can be calculated directly. 6 depends on the atomic forces, but can be estimated 
with sufficient accuracy, since the term containing it is very small. The values are 


y=1.075 x10", 6=—0.56 x 10-™, 
In reality £ is determined from the following formula: 


AM Qo? Leet De nah 
Wie+E @ ((P.- T)+7 (ag) +0 (7e—7)| (11) 


(h, k, l are the Miller indices of the crystal plane), where 


(=) : /() = g2AM 
I }7,-86/ \ I ) 7,=290 ; 


The measurements are made with Mo K, radiation (A = 0.71 A) and the angles involved 
are Up = — Uy, = 1°25, wy = — w, = 0°75. 

To account for the general background scattering James and Brindley ro- 
tated the ionization chamber 2° from the Bragg position and made measurements 
altering the angles of rotation correspondingly. The values thus obtained were sub- 
tracted from the main results. In our calculations we have to make the correspond- 
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cr 


— ———— er 


ileal 


sin? ee 


Ocorr 4 | a 
Sao | a So ee oO. 


T= 290°| T= 86° 


0.178 0.036 0.142 0.037 0.010 0.358 0.098 
0.148 0.035 0.113 0.097 0.026 0.405 0.112 
0.138 | 0.038 0.100 0.129 0.034, 0.398 0.110 
0.129- | 0.036 0.093 0.179 0.048 0.406 0.115 
0.117 0.036 0.081 0.282 0.075 0.391 0.114 
0.115 0.037 0.078 0.298 0.079 0.380 0.111 


Br is . M.v. | 0.390 | 0.110 


-- Table 3 
4 AM 
é Spectr. heir AM Correction AC cone oe 
e+e +P 
i 400 16 0.107 0.013 0.120 0.00750 
600 36 0.242 0.034 0.276 0.00767 
| 444 48 0.344 0.044 0.388 0.00808 
4 800 64 0.450 0.059 0.509 0.00795 
oe 10.00 100 0.714 0.088 0.802 . 0.00802 
; 666 108 0.786 0.093 0.879 0.00814 
a M.V. 0.00789 
a 
“a 
James and Brindley obtain 
; —21 sin” ee 2 
4 B=7.25x10 corresponding to M=1.688 —.,— 4 


After correction of their AM values our result is. 
sin” Up 


je 


B=8.11x10-"! corresponding to M = 1.885 ae 
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. ' : : “1 : = = i 3 a a i = Re a 
| Beikce tiadion: . 7d Wan tawe lft galvtdler br ath ed? sad eee 


‘ < a a — a . -“ 
aa H i$ tea PR? Cee em Te Ae eee 1H ung weit bau esoiohe ae 
pes as At fb, yet Tema fee he atin apie 
Sen Ae Tin tts} a —— 4 E 


ay %. pin . 


The measurements are, however, not very good. Indeed if we use the theo Sai 
t SIMS \aaaee 
values of the atomic scattering factors Fy+and F.- and compute M (=) for 


each temperature, we do not» get the desirable constancy. The constancy of 
AM /(h? +k? +7) is, however, satisfactory, so the value of AM and thus the ex. 
perimental value of M should be better than each of the single measurements.} 


Renninger’s measurements [2] ss 


Renninger has made measurements on rock salt at room temperature and the 
temperature of liquid hydrogen. However, the value of the low temperature seems 
to be very uncertain, so we have only considered the measurements at room temper- 
ature. Renninger has registered the reflection curves and made the integration of 
them. Using the theoretical values of the atomic scattering factors Fy,+ and Fo- 
he calculated the separate temperature factors My,+ and M,- and also determined 
the best common M-value for sodium and chlorine. He also made the reverse calcula- 
tion using the values of Fy,+ and F.-. Plotting these values against sin v)/A 
he obtained curves which agree very well with the theoretical F-curves. Indeed, at 
room temperature they coincide. However, in calculating the characteristic temper- 
ature ©, corresponding to his common M-value, from the well-known formula 


a 
“| sin” vy 


REO +7} 


in?e.\—2 

1 We found that M ("5") at room temperature was decreasing from 2.56 (400) to 
1.75 (666). It was thus constantly higher than the value 1.69 obtained by James and Brindley 
from AM|(h® +k? +1*), At 86° K the value goes from 1.40 down to 0.61 for the same reflections. 
The discrepancy of the single measurements from the expected values can be explained by the 
presence of an absorbing surface layer on the crystal, the possibility of which was pointed 


out by James and Brindley. The lowest reflections must give the largest absorption and in 
calculating AM for each reflection the error is eliminated. 
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tance none is the integral of ihe curve. Thus he aie aie. is more 
uncertain than the value of c. Even in this case the corrections due to the i e 

{ a cancel each other when we take the difference o —o’. one results are given 
Ta 4. 


Table 4 


ae ta) 


| sin? V/A? 


0.0947 0.321 


0.1263 0.306 _ 0.202 0.104 0.024 ; 0.188 
0.2841 0.269 0.166 0.103 0.077 0.261 ~ 
0.3473 0.261 0.157 0.104 0.103 0.282 
0.3789 0.257 0.153 0.104 0.117 0.292 
800 0.5051 0.246 0.143 0.103 0.171 0.313 
644 0.5367 0.244 0.142 0.102 0.183 0.313 
| 822 0.5683 0.242 0.140 0.102 0.197 0.316 
; 10.00 0.7983 0.233 0.129 0.104 0.302 0.334 
_ 666 0.8524 0.231 0.125 0.106 0.337 0.341 
7 = — 
The mean- value: of ‘the eight last results... . - 2. 1.9. - + 2 eee 0.307 
; , 
At room temperature Renninger has obtained 
SUN? Vy ae ; ° 
M =1.15 re A corresponding to ®=319 K. 


Correction for the influence of the diffuse scattering gives 


ait Uy > 4 @ —302°K 
M =1.30 22 A corresponding to © = 30 : 
From measurements of the specific heats we have @=281°K 


and calculations from the elastic constants [9] give ©=303° K. 


The value obtained from specific heats is, however, a mean value over a large 
temperature range, and the value © = 303° K is based upon a Debye frequency distri- 
bution. Thus the perfect agreement between the latter and the value obtained from 
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Renninger’s measurements after correction for the diffuse scattering migh 


ental, but anyhow it is apparent that the high value of 319°K 
by considering the diffuse scattering.’ 


’ Discussion of the results 


The approximation introduced by taking a mean value of the part of 7, that 
depends on the direction of the vector g, means that we replace the real surfaces 
of equal diffusion by spheres. When we integrate over the whole parallelepiped around 
the reciprocal lattice point the errors must cancel each other more or less, but when 
the direction of g varies within a smaller solid angle, as is the case in calculating o’, — 
a larger error may be introduced in the value of the integral. Renninger himself — 
says that the area on his curves corresponding to.A BCD in Fig. 3 amounts to up to 
15 % of the reflection curve, that is at the most 26 % of our calculated Bragg reflection. 
According to our calculations the correction o’ corresponds to up to 40% of the real 
integrated Bragg reflection. : 

The measurements are made without taking account of the azimuth of the crystal. — 
Thus, in Fig. 2 the reciprocal lattice may be turned around OM and the parallelepiped — 
will occupy different regions in the reciprocal lattice. However, Renninger has ob-— 
served that the occasional adjustment of the crystal influences the results of the — 
measurements by about 5%. These facts support the assumption that using the 
mean value does not give any large error in the value of the integrated diffuse — 
reflection, for if that value were very sensitive to the detailed shape of the surfaces 
of equal diffusion within the parallelepiped, the different adjustments of the crystal 
would give corresponding variations of the experimental results. : 

The above calculations from experimental results with allowance for diffuse 
scattering cannot pretend to give the M-values with any great accuracy, because 
the earlier measurements are not very good, and in correcting Renninger’s measure- 
_ments it is difficult to determine how much of the total diffuse scattering influences 

his results. It is evident, however, that the influence of the diffuse scattering is very 
large. Renninger’s high value of the characteristic temperature of rock-salt must 
be due to that circumstance. The zach that the effect of the diffuse scattering can be 
written as a factor exp( const a) explains the good agreement between Ren- 


ninger’s curves of the experimental F'-factors and the curves of the theoretical values. 


IT am most indebted to the head of this institute, Professor Ivar Waller, for suggesting this 
investigation and for much valuable advice during the work. 


Institute of Theoretical Physics, University of Uppsala, Sweden. 


Sisters recent article [10] (which appeared when this paper was going to press) Blackman points 
out that Renninger’s high value of the characteristic temperature of rock-salt cannot be explained 
from the theory of Bragg scattering. 
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